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Q\ ■ We report on a comparative study of the profile of the CKM unitarity triangle, and the resulting 
CP asymmetries in B decays, in the standard model and in several variants of the minimal 
SpH. supersymmetric standard model (MSSM), characterized by a single phase in the quark flavour 
Ph| mixing matrix. The supersymmetric contributions to the mass differences AM^, AM S and to 
the CP-violating quantity |e| are, to an excellent approximation, equal to each other in these 
theories, allowing for a particularly simple way of implementing the resulting constraints on 
the elements of Vqkm from the present knowledge of these quantities. Incorporating the next- 
to-leading-order corrections and applying the current direct and indirect constraints on the 
supersymmetric parameters, we find that the predicted ranges of sin 2/3 in the standard model 
and in MSSM models are very similar. However, precise measurements at 5-factories and 
hadron machines may be able to distinguish these theories in terms of the other two CP- violating 
phases a and 7. This is illustrated for some representative values of the supersymmetric 
contributions in AM^, AM S and lei. 




Figure 1: The unitarity triangle. The angles a, (3 and 7 can be measured via CP violation in 
the B system. 



1 Introduction 



Within the standard model (SM), CP violation is due to the presence of a nonzero complex 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix V [0. A particularly 
useful parametrization of the CKM matrix, due to Wolfenstein |2| , follows from the observation 
that the elements of this matrix exhibit a hierarchy in terms of A, the Cabibbo angle. In this 
parametrization the CKM matrix can be written approximately as 



V 



I l-\X 2 




-X{l+tA 2 X^) 1-|A 2 AX 2 J . (1) 
\AX 3 (1- p-ir]) 



The allowed region in p-r\ space can be elegantly displayed using the so-called unitarity triangle 
(UT). The unitarity of the CKM matrix leads to the following relation: 

v ud v: b + v cd v; b + v td v* = . (2) 

Using the form of the CKM matrix in Eq. (P, this can be recast as 

^+^ = 1, (3) 



\V cb \V a 



which is a triangle relation in the complex plane (i.e. p-rj space), illustrated in Fig. [I]. Thus, 
allowed values of p and rj translate into allowed shapes of the unitarity triangle. 

Constraints on p and rj come from a variety of sources. Of the quantities shown in Fig. 
\V c b\ and \V u b\ can be extracted from semileptonic B decays, while \V t d\ is probed in B®-B® 
mixing. The interior CP- violating angles a, f3 and 7 can be measured through CP asymmetries 
in B decays [[|. Additional constraints come from CP violation in the kaon system (|e|), as 
well as B®-B® mixing. 

In light of the fact that the 5-factory era is almost upon us, one of the purposes of this paper 
is to update the profile of the unitarity triangle within the SM using the latest experimental 
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data. This analysis is done at next-to-leading-order (NLO) precision, taking into account the 
state-of-the-art calculations of the hadronic matrix elements from lattice QCD and available 
data. This therefore provides a theoretically-robust overview of the SM expectations for the 
allowed values of the CP-violating phases, as well as their correlations. 

These CP phases are expected to be measured in the very new future. If Nature is kind, the 
unitarity triangle, as constructed from direct measurements of a, f3 and 7, will be inconsistent 
with that obtained from independent measurements of the sides. If this were to happen, it would 
be clear evidence for the presence of physics beyond the SM, and would be most exciting. 

One type of new physics which has been extensively studied is supersymmetry (SUSY). 
There are a number of hints suggesting that SUSY might indeed be around the corner. One 
example is gauge-coupling unification: a supersymmetric grand unified theory does better than 
its non-supersymmetric counterpart. (How compelling these hints are depends on one's point 
of view.) Partly motivated by this success, but mostly on theoretical grounds, a great deal of 
effort has gone into a systematic study of the pattern of flavour violation in SUSY, where a 
number of flavour- changing neutral-current processes in K and B decays have been studied. 
In this paper we investigate the profile of the unitarity triangle in supersymmetric models. In 
particular, we explore the extent to which SUSY can be discovered through measurements of 
the sides and angles of the unitarity triangle. 

If new physics (of any type) is present, the principal way in which it can enter is via new 
contributions, possibly with new phases, to K°-K°, B®-B® and B®-B® mixing Q. The decay 
amplitudes, being dominated by virtual W exchange, remain essentially unaffected by new 
physics. Thus, even in the presence of new physics, the measured values of |V^| and \V u b\ 
correspond to their true SM values, so that two sides of the UT are unaffected (see Fig. [I]). 
However, the third side, which depends on \V t d\, will in general be affected by new physics. 
Furthermore, the measurements of |e| and B®-B® mixing, which provide additional constraints 
on the UT, will also be affected. Therefore, if new physics is present, the allowed region of the 
UT, as obtained from current experimental data, may not correspond to the true (SM) allowed 
region. 

The CP angles a and (3 are expected to be measured via CP-violating asymmetries in 
B®(t) — > 7r + 7r~ and B®(t) — > J/^/K s , respectively [[|. If there is new physics in B®-B® mixing, 
with new phases, then these measurements will be affected. On the other hand, if there are 
no new phases, then the measurements will probe the true SM values. The third angle 7 can 
be measured in a variety of ways. If it is obtained from a CP asymmetry involving a neutral 
B meson (B® or B®), then it may be affected by new physics. However, if it is measured via 
charged B decays, then the measured value will correspond to the true (SM) value. 

The most general SUSY models allow for the presence of new phases in the couplings of 
supersymmetric and ordinary particles. In such models the new phases are essentially uncon- 
strained, so that the measured CP phases can be greatly shifted from their SM values. If this is 
the case, then the new physics will be relatively easy to find: the unitarity triangle constructed 
from the measurements of a, (3 and 7 will be considerably different from that constructed from 
measurements of the sides. However, precisely because these new phases are unconstrained, 
these models lack predictivity. All that one can say is that there are regions of parameter 
space in which large effects are possible. (Indeed, this is true for any model of new physics 
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which allows for new phases. Examples include models with four generations, Z-mediated 
flavour-changing neutral currents, flavour- changing neutral scalars, etc. ||.) 

In part because of this, most of the theoretical attention has focussed on SUSY models 
which, though less general, have considerably more predictive power. The most theoretically- 
developed model is the minimal supersymmetric standard model (MSSM). Although this is 
often referred to as a single model, in point of fact there are several variants of the MSSM. 
Among these is the scenario of minimal supersymmetric flavour violation ||, which involves, 
in addition to the SM degrees of freedom, charged Higgs bosons, a light stop (assumed right- 
handed) and a light chargino, with all other degrees of freedom assumed heavy and hence 
effectively integrated out. This scenario can be embedded in supergravity (SUGRA) models 
with gauge-mediated supersymmetry breaking, in which the first two squark generations and 
the gluinos are assumed heavy. Regardless of which variant is used, the key point for our 
purposes is that there are no new phases in the couplings - although there are many new 
contributions to meson mixing, all are proportional to the same combination of CKM matrix 
elements as found in the SM. As explained above, in this class of models measurements of the 
CP phases will yield the true SM values for these quantities. However, measurements of meson 
mixing will be affected by the presence of this new physics. 

In this minimal SUSY scenario, NLO calculations for AMd, AM S and |e| can be extracted 
from the work of Krauss and Soff |7]]. Also, NLO corrections to the decay B — > Xgj have 
been worked out by Ciuchini et al. ||. We make use of this work and present the profile 
of the unitarity triangle and CP-violating phases in this scenario, at the NLO precision. A 
particularly nice feature is that the SUSY contributions to K°-K°, B®-B® and B®-B® mixing 
have the same form. Thus, as far as the unitarity triangle is concerned, the various SUSY 
models can be distinguished by a single parameter, /. This simplifies the analysis considerably. 

We note that bits and pieces of such an analysis are already present in the literature. 
However, a theoretically-consistent analysis of the CKM unitarity triangle and the CP-violating 
phases, taking into account all constraints, has not been performed yet, to the best of our 
knowledge. 

With this analysis, one can explore the extent to which the presence of minimal supersym- 
metry can be discovered through the precision measurements of the unitarity triangle which 
will be undertaken by experiments at the 5-factories and hadron colliders. As we will see, the 
profiles of the unitarity triangle in the SM and in MSSM models are similar. However, precise 
measurements may be able to distinguish them. 

In Section 2, we discuss the profile of the unitarity triangle within the SM. We describe the 
input data used in the fits and present the allowed region in p-r\ space, as well as the presently- 
allowed ranges for the CP angles a, (3 and 7. We turn to supersymmetric models in Section 3. 
We review several variants of the MSSM, in which the new phases are essentially zero. We also 
discuss the NLO corrections in such models and show that the SUSY contributions to K°-K°, 
B®-B® and B®-B® mixing are of the same form and can be characterized by a single parameter 
/. We compare the profile of the unitarity triangle in SUSY models, for various values of /, 
with that of the SM. We conclude in Section 4. 
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2 Unitarity Triangle: SM Profile 



2.1 Input Data 

The CKM matrix as parametrized in Eq. ([!]) depends on four parameters: A, A, p and rj. We 
summarize below the experimental and theoretical data which constrain these CKM parameters. 

• \Vus\: We recall that |V^ S | has been extracted with good accuracy from K — > nev and 
hyperon decays || to be 

|K* I = A = 0.2196 ± 0.0023 . (4) 
In our fits, we ignore the small error on A. 

• |Kfe| : The determination of \V c b\ from inclusive and exclusive B decays has been studied 
in a number of papers ||9], |It], [TT], [T2], [T3|]. The number used here is taken from the Particle 
Data Group compilation [§]: 

iKtl = 0.0395 ±0.0017, (5) 

yielding 

A = 0.819 ± 0.035 . (6) 

• |Kfe/Kb| : The knowledge of the CKM matrix element ratio IK&/K&I is based on the anal- 
ysis of the end-point lepton energy spectrum in semileptonic decays B — > X u tvi and the 
measurement of the exclusive semileptonic decays B — > (jr,p)£i/£. Present measurements 
in both the inclusive and exclusive modes are compatible with 



V, 



ub 



V, 



cb 



0.093 ±0.014. (7) 



This gives 



y + v 2 = Q.423 ± 0.064 . (8) 

\e\,B K : The experimental value of |e| is ||: 

|e| = (2.280 ±0.013) x 10~ 3 . (9) 

In the standard model, |e| is essentially proportional to the imaginary part of the box 
diagram for K°-K° mixing and is given by [|TJ]] 

|e| = 6V2^AM K K ( A V ^c{Ws(yc,lfc) -Vcc} 

+ fj tt y t f2(yt)A 2 X\l- P )), (10) 
where t/i = mf/Myy, and the functions / 2 and / 3 are: 

19 1 3 1 3 x 2 lnx 



4 4(l-x) 2(l-x) 2 2(l-x] 
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(The above form for fs(x,y) is an approximation, obtained in the limit i < y. For the 
exact expression, see Ref. |15| .) Here, the fji are QCD correction factors, calculated at 
next-to-leading order in Refs. |L6| (r/ cc ), [T7j (Vtt) and ]Tj| (fj ct ). The theoretical uncer- 
tainty in the expression for |e| is in the renormalization-scale independent parameter Bk, 
which represents our ignorance of the hadronic matrix element (K° 1(^7^(1 — 75) s) \K°). 
Recent calculations of Bk using lattice QCD methods are summarized in Refs. |19|, j20fl , 
yielding 

B K = 0.94 ±0.15. (12) 



AM d , f B< B Bd : The present world average for AM d is plU 

AM d = 0.471 ± 0.016 (ps)- 1 . (13) 

The mass difference AM d is calculated from the B^-B® box diagram. Unlike the kaon 
system, where the contributions of both the c- and the t-quarks in the loop are important, 
this diagram is dominated by t-quark exchange: 

AM d = ^M 2 W M B (f Bd B Bd ) f, B y t f 2 (y t )\V td V tb \ 2 , (14) 

where, using Eq. ([[]), |V^Vt b | 2 = A 2 A 6 [(1 — p) 2 + rf]. Here, fj B is the QCD correction. In 
the fits presented in Ref. |22| we used the value fj B = 0.55, calculated in the MS scheme, 
following Ref. ||17|| . Consistency requires that the top quark mass be rescaled from its 
pole (mass) value of m t = 175 ± 5 GeV to the value rnl{m t {pole)) = 165 ± 5 GeV in the 
MS scheme. We shall ignore the slight dependence of f] B on rnl(m t (pole)) in the range 
given here. 

For the B system, the hadronic uncertainty is given by f B B Bd , analogous to Bk in the 
kaon system, except that in this case f Bd has not been measured. Present estimates of 



this quantity are summarized in Refs. |[T9|, |20| , yielding f Bd y B Bd = (190 ±23) MeV in the 
quenched approximation. The effect of unquenching is not yet understood completely. 
Taking the MILC collaboration estimates of unquenching would increase the central value 



of f Bd JB Bd by 21 MeV |p3| . In our fits, we have taken 



/£ d V^ = 215±40MeV > ( 15 ) 

which is a fairly conservative estimate of the present theoretical error on this quantity. 

AM S , f Bs B Bs : Mixing in the B° s -B® system is quite similar to that in the B d -B^ system. 
The B Q S -B® box diagram is again dominated by t-quark exchange, and the mass difference 
between the mass eigenstates AM S is given by a formula analogous to that of Eq. (|14|) : 

AM S = ^ 2 M 2 w M Bs (f 2 Bs B Bs ) f, Bs y t f 2 (y t )\V*V tb \ 2 . (16) 

Using the fact that \V c b\ = \V ts \ (Eq. (|1|)), it is clear that one of the sides of the unitarity 
triangle, (V^/AV^I, can be obtained from the ratio of AM d and AM S , 



AM S _ Vb s M Bs {f 2 Ba B l 



AM, fj Bd M Bd (f 2 B Bc 



V ts \ 2 



Vt 



td 



(17) 
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The only real uncertainty in this quantity is the ratio of hadronic matrix elements 

fB s BB s / fl d BB d - It is now widely accepted that the ratio £ s = ( }b s \[b~b~ s ) / (/s d V B Bd ) 
is probably the most reliable of the lattice-QCD estimates in B physics, and the present 
estimate is [J9|, fZOfl : 

£ s = 1.14 ±0.06. (18) 



The present lower bound on AM S is: AM S > 12.4 (ps) (at 95% C.L.) (TJ. This bound 



has been established using the so-called "amplitude method" |24]]. This method is also 
the best procedure for including information about B° s -B® mixing in the fit, and works as 
follows. Given a meson which at t = was pure B° s) the probability for it to be detected 
as a B° s at time t is 

P = ^-e~ t/T (1 + cos AM s t) , (19) 
It 

while the probability for it to be detected as a B® is 

P = ^e~ t/T (1 - cos AM s t) . (20) 
2r 

In the amplitude method, one introduces the amplitude A and writes the oscillations 
terms as (1 ± ^4 cos AM s t). One then measures the value of A, along with its error a a 
assuming various values of AM S . For a given value of AM S , if A is compatible with 0, 
one concludes that there is no visible oscillation at this frequency; if A is compatible with 
1, one concludes that an oscillation was observed at this frequency. 

The experimental data consists of measured values of A and for various values of AM S . 
To include this data in the fit, for each set of free parameters (A, p, 77, £ s ) we calculate 
the value of AM S and find the corresponding experimental values of A and a^. Since a 
nonzero value of AM S implies that there is B° s -B® mixing, theoretically one should have 
A = 1. For this set of parameters we therefore add to the global \ 2 a factor 

A- 1\ 2 

±-±) . (21) 

(There is also a similar factor which takes into account the deviation of £ s from its central 
value (Eq. ©).) 



The data used in our fits are summarized in Table |]. The quantities with the largest errors 

are f] cc (28%), B K (16%), \V ub /V cb \ (15%) and f Bi 4 B ~B i ( 19 %)- Of these, the latter three are 
extremely important in defining the allowed p-rj region (the large error on r) cc does not affect 

the fit very much). The errors on two of these quantities — B K and fs d \J B Bd — are purely 
theoretical in origin, and the error on |V^b/V^,| has a significant theoretical component (model 
dependence). Thus, the present uncertainty in the shape of the unitarity triangle is due in large 
part to theoretical errors. Reducing these errors will be quite important in getting a precise 
profile of the unitarity triangle and the CP-violating phases. 

There are two other measurements which should be mentioned here. 
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Parameter 


Value 


A 


0.2196 




v cb \ 


0.0395 ±0.0017 




v ub /v cb \ 


0.093 ±0.014 




e 


(2.280 ±0.013) x lO" 3 


AM d 


(0.471 ±0.016) (ps)- 1 


A M s 


> 12.4 (ps)- 1 


rrTi(rnt(pole)) 


(165 ±5) GeV 


m^{m c (pole)) 


1.25 ±0.05 GeV 


Vb 


0.55 


Vcc 


1.38 ±0.53 


Vet 


0.47 ±0.04 


Vtt 


0.57 


Br 


0.94 ±0.15 


fB d \jB Bd 


215 ±40 MeV 


6 


1.14 ±0.06 



Table 1: Data used in the CKM fits. Values of the hadronic quantities Bk, fB d \BB d , and 
£ s are taken from the lattice QCD results [19], BQ]. The remaining theoretical numbers are 



discussed in the text. The value for AM^ and the 95% CL. lower limit on AM S are taken from 
the LEP Electroweak group 



12 1 . All other experimental numbers are taken from the Particle 



Data Group 
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First, the KTEV collaboration [25| has recently reported a measurement of direct CP vio- 
lation in the K sector through the ratio e'/e, with 



Re(eVe) = (28.0 ± 3.0(stat) ± 2.6(syst) ± 1.0(MC stat)) x 10" 



(22) 



in agreement with the earlier measurement by the CERN experiment NA31 ||26|| , which reported 
a value of (23 ± 6.5) x 10~ 4 for the same quantity. The present world average is Re(e'/e) = 
(21.8 ±3.0) x 10 -4 . This combined result excludes the superweak model [27 by more than 7a. 



A great deal of theoretical effort has gone into calculating this quantity at next-to-leading 
order accuracy in the SM [28], |29|, ROf. The result of this calculation has been summarized in a 



succint form by Buras and Silvestrini [31]: 



Re (e'/e) = ImA t 



1.35 + iL (l.l|ri 8) |Sj 1/2) + (1.0 



0.67|r^|)^ 3/2) ) 



(23) 



Here \t = VtdV t * = A 2 X 4 7] and r z represents the short-distance contribution, which at the 
NLO precision is estimated to lie in the range 6.5 < 



,(8) I 



< 8.5 28 



(1/2) 



B, 



(1/2), 



and B 



(3/2) 



B 



(3/2), 



8 V"V 

AI = 3/2 operators 6 and 8 , respectively, calculated at the scale \i 
and the l/N c expansion |33) yield: 



J. The quantities 
are the matrix elements of the AI = 1/2 and 

m c . Lattice-QCD [[32 



0.8 < B, 



(1/2) 
6 



< 1.3, 



0.6 < B 



(3/2) 



< 1.0 . 



Finally, the quantity R s in Eq. ( [23] ) is defined as: 

( 150 MeV 

Re = 



m s (m c ) + m d (m c ) 



(24) 



(25) 



essentially reflecting the s-quark mass dependence. The present uncertainty on the CKM 
matrix element is ±23%, which is already substantial. However, the theoretical uncertainties 
related to the other quantities discussed above are considerably larger. For example, the ranges 
e'/e = (5.3 ± 3.8) x 10~ 4 and e'/e = (8.5 ± 5.9) x 10~ 4 , assuming m s (m c ) = 150 ± 20 MeV 
and m s (m c ) = 125 ± 20 MeV, respectively, have been quoted as the best representation of the 
status of e'/e in the SM |3~4j| . These estimates are somewhat on the lower side compared to the 
data but not inconsistent. (For some recent speculations on new physics, see Refs. [|5|, |36fl.) 

Thus, whereas e'/e represents a landmark measurement, removing the superweak model of 
Wolfenstein and its kith and kin from further consideration, its impact on CKM phenomenology, 
particularly in constraining the CKM parameters, is marginal. Probably the best use of this 
measurement is to constrain the s-quark mass, which at present has considerable uncertainty. 
For this reason we do not include the measurement of e'/e in the CKM fits presented here. 

Second, the CDF collaboration has recently made a measurement of sin 2(3 |37| |. In the 
Wolfenstein parametrization, —f3 is the phase of the CKM matrix element V t d- From Eq. (|T|) 
one can readily find that 

- 0^ • < 26 > 
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Figure 2: Allowed region in p-rj space in the SM, from a fit to the ten parameters discussed 
in the text and given in Table [I]. The limit on AM S is included using the amplitde method of 



Ref. [p4| . The theoretical errors on f Bd y BB d , Bk and £<, are treated as Gaussian. The solid line 
represents the region with % 2 = Xmin + 6 corresponding to the 95% C.L. region. The triangle 
shows the best fit. 

Thus, a measurement of sin 2/3 would put a strong contraint on the parameters p and rj. How- 
ever, the CDF measurement gives P?| 

sin2/? = 0.79_tg;£ , (27) 

or sin 2(3 > at 93% C.L. As we will see in the next section, this constraint is quite weak - the 
other measurements already constrain 0.52 < sin 2/3 < 0.94 at the 95% C.L. in the SM. (The 
CKM fits reported in Refs. |12|, [13|, [38| yield similar ranges.) In light of this, and given that 



it is not clear how to combine the above measurement (which allows for unphysical values of 
sin 2(3) with the other data, we have not included this measurement in our fits. 



2.2 SM Fits 



In order to find the allowed region in p-rj space, i.e. the allowed shapes of the unitarity triangle, 
the computer program MINUIT is used to fit the parameters to the constraints described above. 

In the fit, we allow ten parameters to vary: p, rj, A, m t , m c , r] cc , r] ct , js d \J Bs d , Bk, and £ s . 
The AM 8 constraint is included using the amplitude method. 

The allowed (95% C.L.) p-rj region is shown in Fig. ||. The best fit has (p, rj) = (0.20, 0.37). 

There is an alternative way to include the AM S constraint, one which we have used in the 
past p2 |. In this case the constraint is excluded from the fit, but rather one cuts away a region 



of p-rj space by superimposing a line corresponding to a particular value of £ s . Using the 95% 
C.L. limit on AM S and allowing AM^ to fluctuate upward by la from its central value, one 
obtains AM S / AM^ > 25.5. This yields the following constraint on p and rj: 



A 



P? + v 2 



( AM d 
V AM, 



25.5 



(1.14±0.06) s 



(28) 
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Figure 3: Allowed region in p-rj space in the SM. The AM S limit is not included in the fit, but 
rather imposed by cutting away a region of the space. The disallowed region is shown for three 
choices of the parameter £ s : 1.08 (dotted line), 1.14 (dashed line), 1.20 (solid line). In all cases, 
the region to the left of the curve is ruled out. 

We take three candidate values for £ s : 1.08, 1.14, 1.20. The results are shown in Fig. |3|. The 
best fit has (p,r?) = (0.19,0.37). 

A comparison of Figs. ^ and |3] reveals that the allowed regions are similar, though not 
identical, for the two methods of including the AM, constraint. Henceforth, all our results will 
be presented using the fit which includes this constraint via the amplitude method. 

The CP angles a, (3 and 7 can be measured in CP- violating rate asymmetries in B decays 
0. For example, the asymmetries in B®(t) — > 7r + 7r~ and B^(t) — > J/^>K S probe sin 2a and 
sin 2/5, respectively. The angle 7 can be extracted from B ± — > DK ± |3{| or B° s {t) — > DfK T 
f40f . The function in this case is sin 2 7. 

From Fig. [I], it is clear that these angles can be expressed in terms of p and r\. Thus, different 
shapes of the unitarity triangle are equivalent to different values of the CP angles. Referring to 
Fig. ||, we note that the preferred (central) values of these angles are (a, /3, 7) = (93°, 25°, 62°). 
The allowed ranges at 95% C.L. are 

65° < a < 123° 
16° < < 35° 

36° < 7 < 97° (29) 

or, equivalently, 

- 0.91 < sin 2a < 0.77 
0.52 < sin 2(3 < 0.94 

0.35 < sin 2 7 < 1.00 . (30) 

Of course, the values of a, (3 and 7 are correlated, i.e. they are not all allowed simultaneously. 
After all, the sum of these angles must equal 180°. We illustrate these correlations in Figs. |] 
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= 215±40MeV, B =0.94+0.15, f = 



0.9 
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= 215±40MeV, B =0.94+0.15, f = 0.75 
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Figure 4: Allowed 95% CL. region of the CP- violating quantities sin 2a and sin 2f3, from a fit 
to the data given in Table [[]. The upper left plot (/ = 0) corresponds to the SM, while the 
other plots (/ = 0.2, 0.4, 0.75) correspond to various SUSY models. 

and |5|. The plots with / = correspond to the SM. (Nonzero values of / correspond to 
supersymmetric models, which will be discussed in detail in the next section.) Fig. |] shows 
the allowed region in sin 2a-sin 2(3 space allowed by the data. And Fig. |^ shows the allowed 
(correlated) values of the CP angles a and 7. This correlation is roughly linear, due to the 
relatively small allowed range of (3 (Eq. (p9|)). 



3 Unitarity Triangle: A SUSY Profile 

In this section we examine the profile of the unitarity triangle in supersymmetric (SUSY) 
theories. As we will see below, the most general models contain a number of unconstrained 
phases and so are not sufficiently predictive to perform such an analysis. However, there is a 
class of SUSY models in which these phases are constrained to be approximately zero, which 
greatly increases the predictivity. Not surprisingly, the calculations of the SUSY contributions 
to measured quantities are also more advanced in these models, having been performed at the 
next-to-leading-order (NLO) level. In the following subsections, we discuss aspects of more 
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Figure 5: Allowed 95% C.L. region of the CP-violating quantities a and 7, from a fit to the 
data given in Table |l[ The upper left plot (/ = 0) corresponds to the SM, while the other plots 
(/ = 0.2, 0.4, 0.75) correspond to various SUSY models. 
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general SUSY theories, as well as the details of that class of theories whose effects on the 
unitarity triangle can be directly analyzed. 



3.1 Flavour Violation in SUSY Models - Overview 

We begin with a brief review of flavour violation in the minimal supersymmetric standard model 
(MSSM). 

The low energy effective theory in the MSSM can be specified in terms of the chiral su- 
perfields for the three generations of quarks (Qi, U?, and Df) and leptons (Lj and Ef), chiral 
superfields for two Higgs doublets {H\ and H 2 ), and vector superfields for the gauge group 
SU(3) C x S77(2)/ x U{l) Y [0. The superpotential is given by 



W MSSM = fZQiDjH! + fjQiUjHz + ftEiLjHi + ^H X H 2 . (31) 

The indices i,j = 1,2,3 are generation indices and f^, f^J, fl are Yukawa coupling matrices 
in the generation space. A general form of the soft SUSY-breaking term is given by 

- £ soft = (m%))q^ + (m 2 D ))d^ + (m^u* + {m 2 E ))e^ + {m\)%P 
+ l\\h\h x + Alh\h 2 - (Biih 1 h 2 + h.c.) 
+ (A^qidjhi + A l jjqiUjh 2 + A l [edjhi + h.c.) 

+ (^BB + ^WW + ^GG + h.c), (32) 

where qi, iti, di, £i, e«, hi and h 2 are scalar components of the superfields Qi, Ui, Di, Li, E iy 
Hi and H 2 , respectively, and B, W and G are the U(l), SU(2) and SU(3) gauge fermions, 
respectively. 

In fact, there are many variants of the MSSM. All have the same particle content, but the 
mass hierarchies and supersymmetry-breaking parameters are different. As a result, they lead 
to different low-energy predictions. 

For example, in supergravity (SUGRA) models, the SUSY-breaking parameters in the 
MSSM are assumed to have a simple structure at the Planck scale ||41|| . Neglecting the dif- 
ference between the Planck and the GUT scales, and defining this scale by X, one can express 
them as follows: 





= (ml)) = m 2 5}, 






= K)J = (ml)) = 


ml5), 


A? 






A% 


= f l DX A xm , A l l = 


= fLx A xm 


Mi 


= M 2 = M 3 = M gX 





(33) 

In the minimal SUGRA case, m = A , while in the non-minimal SUGRA case, one takes m 
and A as independent parameters. 
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In general, MSSM models have three physical phases, apart from the QCD vacuum param- 
eter Q QCD which we shall take to be zero. The three phases are: (i) the CKM phase represented 
here by the Wolfenstein parameter rj, (ii) the phase 9 a = ar g(^4), and (iii) the phase 9^ = arg(/x) 
42f . The last two phases are peculiar to SUSY models and their effects must be taken into 



account in a general supersymmetric framework. In particular, the CP-violating asymmetries 
which result from the interference between mixing and decay amplitudes can produce non- 
standard effects. Concentrating here on the AB = 2 amplitudes, two new phases 9d and 9 S 



arise, which can be parametrized as follows [43 



7Js 




(Bd, s \ 


USUSY 


\Bd,s)\ 


{Bd, s 


\ n eff\ 


Bd,s) ) 



(34) 



where r ri SUSY is the effective Hamiltonian including both the SM degrees of freedom and the 
SUSY contributions. Thus, CP-violating asymmetries in B decays would involve not only the 
phases a, (3 and 7, defined previously, but additionally 6^ or 9 S . In other words, the SUSY 
contributions to the real parts of Mi 2 (Bd) and Mi 2 (B s ) are no longer proportional to the CKM 
matrix elements V t dV t l and Vt s V t * b , respectively. If 9a or 9 S were unconstrained, one could not 
make firm predictions about the CP asymmetries in SUSY models. In this case, an analysis of 
the profile of the unitarity triangle in such models would be futile. 

However, the experimental upper limits on the electric dipole moments (EDM) of the neu- 



tron and electron || do provide constraints on the phases 9^ and 9a @- In SUGRA models 
with a priori complex parameters A and fi, the phase 9^ is strongly bounded with 9^ < 0.017T 
|[45|| . The phase 9a can be of 0(1) in the small 9^ region, as far as the EDMs are concerned. 
In both the AS = 2 and AB = 2 transitions, and for low-to- moderate values of tanu [], it 
has been shown that 9a does not change the phase of either the matrix element Mi2(K) 
or of M 12 (B) ||. Hence, in SUGRA models, arg M 12 (B)\ SUGRA = argM 12 (£)| SM = arg(& 2 ), 
where & = V t * d V t b. Likewise, the phase of the SUSY contribution in M\ 2 (K) is aligned with the 
phase of the tt-contribution in Mi 2 (K), given by arg(V^V^*). Thus, in SUGRA models, one can 
effectively set 9^ — and 9 S ~ 0, so that the CP-violating asymmetries give information about 
the SM phases a, (3 and 7. Hence, an analysis of the UT and CP-violating phases a, (3 and 7 
can be carried out in a very similar fashion as in the SM, taking into account the additional 
contributions to M 12 (K) and M 12 (-E>). 

For large-tan v solutions, one has to extend the basis of H e ff(AB = 2) so as to include new 
operators whose contribution is small in the low-tan v limit. The resulting effective Hamiltonian 
is given by 

C 2 M 2 3 

H eff (AB = 2) = -^—w £ cMOi , (35) 

i=l 

where O x = dlj^bld^b^, 2 = d\b%d{b P R and 3 = d\\P R d\b^ and d are the Wilson co- 



efficients |46|, [47j. The coefficients Ci(fi) and C 2 (fx) are real relative to the SM contribution. 



supersymmetric jargon, the quantity tan/3 is used to define the ratio of the two vacuum expectation 
values (vevs) tan/3 = v u /vd, where Vd(v u ) is the vev of the Higgs field which couples exclusively to down- type 
(up-type) quarks and leptons. (See, for example, the review by Haber in Ref. ||). However, in discussing 
flavour physics, the symbol (3 is traditionally reserved for one of the angles of the unitarity triangle. To avoid 
confusion, we will call the ratio of the vevs tanu. 
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However, the chargino contributions to C^(fi) are generally complex relative to the SM contri- 
bution and can generate a new phase shift in the B°-B° mixing amplitude |48], 49j . This effect 



is in fact significant for large tant; j46], since C 3 (/i) is proportional to {m b /m w cos (3) 2 . How 
large this additional phase (9 a and 9 S ) can be depends on how the constraints from EDM are 
imposed. For example, Baek and Ko |39[] find that in the MSSM without imposing the EDM 



constraint, one has 2\9d\ < 6° for a light stop and large tan-u but this phase becomes practically 
zero if the EDM constraints [OT are imposed. This aspect of the analysis in Ref. |49|j , though 
done without invoking the SUGRA model mass relations, should also hold in SUGRA models. 

In a more general supersymmetric scenario, the SUSY effects may lead to additional phases 
and the constraints from EDMs may not be sufficient to effectively bound them. (A recent 
example is the so-called "effective super symmetry" |43| .) As discussed above, such models do 
not make clean predictions for the CP asymmetries, so that one cannot analyse the profile of 
the unitarity triangle in this case. Instead, flavour- changing effects in K and B decays in these 
scenarios will have to be disentangled along the lines suggested in Refs. J|3], [51], |52| . These 
theories may give rise to flavour- violation effects which are uncharacteristic of the SM and 
MSSM. In particular, they may lead to a measurable charge asymmetry in semileptonic decays 
of the B° and W |>2|, which is estimated to be negligible in both the SM and MSSM/SUGRA 



models discussed above. However, some technical aspects in the most general supersymmetric 
theories remain to be worked out. For example, the complete NLO corrections for AM^, 
AM S and |e| have not been calculated in these theories. The bag constants which correct for 
the vacuum insertion approximation in the matrix elements of the operators in H^ =2 have 
been calculated using lattice-QCD techniques [(53j, but the corresponding estimates for the 
bag parameters H^ =2 are not yet available. Likewise, the NLO calculations in the general 
supersymmetric case for the decay B — > X s ^ are not yet at hand. This last ingredient enters 
the estimates of the quantities AM^, AM S and |e| indirectly as the measured branching ratio 
BiB — > X s 7) provides rather stringent constraints on the allowed supersymmetric parameters. 

In view of the foregoing, we shall restrict ourselves to a class of SUSY models in which the 
following features, related to flavour mixing, hold: 



• The squark flavour mixing matrix which diagonalizes the squark mass matrix is approx- 
imately the same as the corresponding quark mixing matrix Vckm, apart from the left- 
right mixing of the top squarks. 

• The first- and second-generation squarks with the same gauge quantum numbers remain 
highly degenerate in masses but the third-generation squarks, especially the top squark, 
can be significantly lighter due to the renormalization effect of the top Yukawa coupling 
constants. 

• The phases 9 d and 9 S are negligible in the entire tant> plane, once the constraints from 
the EDMs of neutron and lepton are consistently imposed. 



These features lead to an enormous simplification in the flavour structure of the SUSY con- 
tributions to flavour- changing processes. In particular, SUSY contributions to the transitions 
b — ► s, b — ► d and s — > d are proportional to the CKM factors, V t bV t * s , V t bV t * d and V ts V* d , re- 
spectively. Similarly, the SUSY contributions to the mass differences Mi 2 (B s ), Mi 2 (-Bd) and 
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M l2 (K) are proportional to the CKM factors (V tb V*) 2 , (V tb V t * d ) 2 and (V ts V t * d ) 2 , respectively. 
These are precisely the same factors which govern the contribution of the top quark in these 
transitions in the standard model. Thus, the supersymmetric contributions can be implemented 
in a straightforward way by adding a (supersymmetric) piece in each of the above mentioned 
amplitudes to the corresponding top quark contribution in the SM. 



3.2 NLO Corrections to AM d , AM S and e in Minimal SUSY Flavour 
Violation 

A number of SUSY models share the features mentioned in the previous subsection, and the 
supersymmetric contributions to the mass differences M\ 2 {B) and Mi 2 (K) have been analyzed 
in a number of papers [[45], |54|, [55], [56], [57j, following the pioneering work of Ref. |58[ . 



Following these papers, AM d can be expressed as: 

AM d = ^M^M B (f 2 Bd B Bd ) f) B [A SM (B) + A H ±{B) + A X± {B) + A- g (B)] , (36) 

where the function Asm(B) can be written by inspection from Eq. (|l4j): 

A SM (B) = y t f 2 (y t )\V t * d V tb \ 2 . (37) 

The expressions for Ajj±(B), A X ±(B) and Ag(B) are obtained from the SUSY box diagrams. 
Here, H ± , xfi ta and dj represent, respectively, the charged Higgs, chargino, stop and down- 
type squarks. The contribution of the intermediate states involving neutralinos is small and 
usually neglected. The expressions for A H ±(B), A X ±(B) and Ag(B) are given explicitly in the 



literature H6|, M, EB 



We shall not be using the measured value of the mass difference AM K due to the uncertain 
contribution of the long-distance contribution. However, |e| is a short-distance dominated 
quantity and in supersymmetric theories can be expressed as follows: 

|e| = 6 ^2 amJ ^ [Im AsM{K) + Im Ah±{K) + Im Ax±(K) + Im MK)] ' (38) 
where, again by inspection with the SM expression for |e| given in Eq. (|T0|), one has 

Im Asm(K) = A 2 X & V {y c {fj ct f 3 (y c , y t ) - fj cc } + fj a y t f 2 (y t )A 2 X\l - p)) (39) 

The expressions for Im A B ±(K), Im A X ±(B) and Im A g (B) can be found in Refs. |^6, ^ . 



For the analysis reported here, we follow the scenario called minimal flavour violation in 
Ref. |J. In this class of supersymmetric theories, apart from the SM degrees of freedom, only 
charged Higgses, charginos and a light stop (assumed to be right-handed) contribute, with all 
other supersymmetric particles integrated out. This scenario is effectively implemented in a 
class of SUGRA models (both minimal and non- minimal) and gauge-mediated models in 



which the first two squark generations are heavy and the contribution from the intermediate 



gluino-squark states is small [[|5|, |54|, [55], |56|, |57 . 



16 



For these models, the next-to-leading-order (NLO) corrections for AM d , AM S and |e| can 
be found in Ref. J7|. Moreover, the branching ratio B{B — > Xg'y) has been calculated in 
Ref. ||. We make use of this information and quantitatively examine the unitarity triangle, 
CP-violating asymmetries and their correlations for this class of supersymmetric theories. The 
phenomenological profiles of the unitarity triangle and CP phases for the SM and this class of 
supersymmetric models can thus be meaningfully compared. Given the high precision on the 
phases a, [3 and 7 expected from experiments at S-factories and hadron colliders, a quantitative 
comparison of this kind could provide a means of discriminating between the SM and this class 
of MSSM's. 

The NLO QCD-corrected effective Hamiltonian for AB = 2 transitions in the minimal 
flavour violation SUSY framework can be expressed as follows [0: 

H eff = ^ 2 (V td V t l)%AB)SO LL , (40) 
where the NLO QCD correction factor rj 2 s(B) is given by 0: 



a s (m w ) (D 



(41) 



AM d = ^(V td V t l)%, s (B)S(f 2 Bd B Bd ) . (42) 



in which n/ is the number of active quark flavours (here n/ = 5), and the quantities Z Hf , 7^ 
and (3^1 are defined below. The operator Oll — 0\ is the one which is present in the SM, 
previously defined in the discussion following Eq. (|35|). The explicit expression for the function 
S can be obtained from Ref. |j8| and for D it is given in Ref. [[/J, where it is derived in the 
NDR (naive dimensional regularization) scheme using MS'-renormalization. 

The Hamiltonian given above for B^S® mixing leads to the mass difference 

Gj 
6tt' z 

The corresponding expression for AM S is obtained by making the appropriate replacements. 
Since the QCD correction factors are identical for AM d and AM S , it follows that the quantities 
AMd an d AM S are enhanced by the same factor in minimal flavour violation supersymmetry, 
as compared to their SM values, but the ratio AM s /AM d in this theory is the same as in the 
SM. 

The NLO QCD-corrected Hamiltonian for AS* = 2 transitions in the minimal flavour viola- 
tion supersymmetric framework has also been obtained in Ref. J7|. From this, the result for e 
can be written as: 

+ UK)SA 2 \\l-p)), (43) 
where the NLO QCD correction factor is 0: 

\a s {m c ) J \ a s {m b ) ) 

1 + ^>(Z 3 - Z0 + ^(Z, - Z.) + =^(§ + zA . (44) 
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Here 

<v(0) 



Z = J^l 1 d (1) (Ah) 

and the quantities entering in Eqs. fl4T|) and (f4"4|) are the coefficients of the well-known beta 
function and anomalous dimensions in QCD: 

1 Af ' 3 



c 



= y ^ - f Nc n f - 2C f n f 



f 

(i) 

l f 2N. 



N c -1 



c 



ni 57 19 Ar 4 
-2H N c + -n f 

N c 3 3 i 



(46) 



with N c = 3 and Cp = 4/3. The ratio 

fj2, s (B)(NLO) = a s (M w ) D 
f) 2 , s (B)(LO) An l S 



+ (47) 



is worked out numerically in Ref. [0 as a function of the supersymmetric parameters (chargino 
mass m X2 , mass of the lighter of the two stops rni R , and the mixing angle in the stop sector). 
This ratio is remarkably stable against variations in the mentioned parameters and is found 
numerically to be about 0.89. Since in the LO approximation the QCD correction factor 
rj 2 ,s(B)(LO) is the same in the SM and SUSY, with 

r) 2 , s (B)(LO) = a s (M w y {0)/2A °} , (48) 

the QCD correction factor f)2 ; s{B)(NLO) entering in the expressions for AM^ and AM S in the 
MSSM is found to be f)2,s(B)(NLO) = 0.51 in the MS'-scheme. This is to be compared with 
the corresponding quantity t)b = 0.55 in the SM. Thus, NLO corrections in AM^ (and AM S ) 
are similar in the SM and MSSM, but not identical. 

The expression for rj2 t s{K)(NLO)/f}2 t s{K)(LO) can be expressed in terms of the ratio 
V2,s{B)(NLO)/f)2,s(B)(LO) given above and the flavour- dependent matching factors Z nf : 

r)2,s{K){LO) r]2,s{B){LO) An An 

where we have used the numerical value fj2,si.B){NLO) /fj 2 ,s{B)(LO) = 0.89 calculated by 
Krauss and Soff 0, along with a s {m c ) = 0.34 and a s (mb) = 0.22. Using the expression for 
the quantity r)2,s{K)(LO), which is given by the prefactor multiplying the square bracket in 
Eq. (p4|), one gets fj 2; s(K)(NLO) = 0.53 in the MS'-scheme. This is to be compared with the 
corresponding QCD correction factor in the SM, fj tt = 0.57, given in Table [l|. Thus the two 
NLO factors are again very similar but not identical. 

Following the above discussion, the SUSY contributions to AM d , AM S and |e| in supersym- 
metric theories are incorporated in our analysis in a simple form: 

AM d = AM d (SM)[l + f d (rn x ±,m iR ,m H ±,tBJiv)}, 
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AM, = AM S (SM)[1 + f s (m x ±,m iR ,m H ±,tanv)}, 

+ fi tt y t f 2 {yt)[l + fe(m x ±,m h ,m H ±,t&nv)]A 2 \ 4 (l - p)). (50) 
The quantities fa, f s and fa can be expressed as 

/„ = /. = H^RA d (S), 
Vb 

h = ^^RaAS), (51) 
Vtt 

where R,A d {S) is defined as 

AM d (SUSY) S 
= AMASM) {W) - ■ <52) 

The functions fa, i — d,s,e are all positive definite, i.e. the supersymmetric contributions add 
constructively to the SM contributions in the entire allowed supersymmetric parameter space. 
We find that the two QCD correction factors appearing in Eq. fl5T] ) are numerically very close 
to one another, with rj2,s(B)/r) B ~ fj2,s(K)/Vtt — 0.93. Thus, to an excellent approximation, 
one has fa = fa = fa = f. 

How big can / be? This quantity is a function of the masses of the top squark, chargino 
and the charged Higgs, m t - fl , m x ± and m^±, respectively, as well as of tan?;. The maximum 
allowed value of / depends on the model (minimal SUGRA, non-minimal SUGRA, MSSM with 
constraints from EDMs, etc.). We have numerically calculated the quantity / by varying the 
SUSY parameters 0, m t - fl , fn X2 , m H ± and tan-u. Using, for the sake of illustration, m x ± = 
mi R = mn± = 100 GeV, m^± = 400 GeV and tan-u = 2, and all other supersymmetric masses 
much heavier, of 0(1) TeV, we find that the quantity / varies in the range: 

0.4 <f < 0.75 for |0| < vr/4 , (53) 

with the maximum value of / being at <fi = 0. These parametric values are allowed by the 
constraints from the NLO analysis of the decay B —>■ X s + 7 reported in Ref. 0, as well as 
from direct searches of the supersymmetric particles 0. The allowed range of / is reduced 
as tan-u increases. Thus, for tant> = 4, one has 0.15 < / < 0.42 for \<p\ < tt/A. Likewise, / 
decreases as m^ R , m X2 and m H ± increase, though the dependence of / on m H ± is rather mild 
due to the compensating effect of the and chargino contributions in the MSSM, as observed 
in Ref. |J. This sets the size of / allowed by the present constraints in the minimal flavour 
violation version of the MSSM. 

If additional constraints on the supersymmetry breaking parameters are imposed, as is the 
case in the minimal and non-minimal versions of the SUGRA models, then the allowed values 
of / will be further restricted. A complete NLO analysis of / would require a monte-carlo 
approach implementing all the experimental and theoretical constraints (such as the SUGRA- 
type mass relations). In particular, the NLO correlation between B(B — ► Xg'j) and / has to 
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be studied in an analogous fashion, as has been done, for example, in Refs. [56, 57] with the 
leading order SUSY effects. 

In this paper we adopt an approximate method to constrain / in SUGRA-type models. We 
take the maximum allowed values of the quantity R,A d (S), defined earlier, from the existing LO 
analysis of the same and obtain / by using Eq. (|5lf) . For the sake of definiteness, we use the 
updated work of Goto et al. |56|, |57J], which is based, among other constraints, on the following: 

• The lightest chargino mass is larger than 91 GeV, and all other charged SUSY particle 
masses are larger than 80 GeV. 

• The gluino and squark masses are bounded from the searches at TEVATRON and LEP 



• Constraints on the supersymmetric parameters taking into account the NLO calculation 
of the decay branching ratio BiB — > X s 7) in the SM and the charged Higgs contribution 
in the MSSM ||60| , and the updated experimental branching ratio from the CLEO and 
ALEPH collaborations, which at 95% C.L. is given by 

2.0 x 10" 4 < B(B -> X s7 ) < 4.5 x 10~ 4 . (54) 

This last constraint plays a rather crucial role in determining the allowed values of R& d (S) 
- and hence of AM d (SUSY)/AM d (SM) and \e\(SUSY)/\e\(SM) — as the magnitude of the 
SUSY contribution in AM d , |e| and B(B — > X s 7) are strongly correlated for the small tan-u 
case [^, |57j. In fact, were it not for the bounds on B[B — > X s ^) given above, much larger 
values of RA d (S) would be allowed. 

From the published results we conclude that typically / can be as large as 0.45 in non- 
minimal SUGRA models for low tant; (typically tanu = 2) |5(|, and approximately half of 
this value in minimal SUGRA models |45|, Relaxing the SUGRA mass constraints, 

admitting complex values of A and \i but incorporating the EDM constraints, and imposing 



the constraints mentioned above, Baek and Ko [[49] find that / could be as large as / = 0.75. 
In all cases, the value of / decreases with increasing tan-u or increasing m^± and rrif R , as noted 
above. 



3.3 SUSY Fits 

For the SUSY fits, we use the same program as for the SM fits, except that the theoretical 
expressions for AM d , AM S and |e| are modified as in Eq. flSUp. We compare the fits for four 
representative values of the SUSY function / — 0, 0.2, 0.4 and 0.75 — which are typical of 
the SM, minimal SUGRA models, non-minimal SUGRA models, and non-SUGRA models with 
EDM constraints, respectively. 

The allowed 95% C.L. regions for these four values of / are all plotted in Fig. || As is clear 
from this figure, there is still a considerable overlap between the / = (SM) and / = 0.75 
regions. However, there are also regions allowed for one value of / which are excluded for 
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Figure 6: Allowed 95% CL. region in p-rj space in the SM and in SUSY models, from a fit to 
the data given in Table |l|. From left to right, the allowed regions correspond to / = (SM, 
solid line), / = 0.2 (long dashed line), / = 0.4 (short dashed line), / = 0.75 (dotted line). 

another value. Thus a sufficiently precise determination of the unitarity triangle might be able 
to exclude certain values of / (including the SM, / = 0). 

From Fig. ^| it is clear that a measurement of the CP angle (3 will not distinguish among the 
various values of /: even with the naked eye it is evident that the allowed range for (5 is roughly 
the same for all models. Rather, it is the measurement of 7 or a which has the potential to rule 
out certain values of /. As / increases, the allowed region moves slightly down and towards 
the right in the p-r\ plane, corresponding to smaller values of 7 (or equivalently, larger values 
of a). We illustrate this in Table 0, where we present the allowed ranges of a, (3 and 7, as 
well as their central values (corresponding to the preferred values of p and 77), for each of the 
four values of /. From this Table, we see that the allowed range of f3 is largely insensitive to 
the model. Conversely, the allowed values of a and 7 do depend somewhat strongly on the 
chosen value of /. Note, however, that one is not guaranteed to be able to distinguish among 
the various models: as mentioned above, there is still significant overlap among all four models. 
Thus, depending on what values of a and 7 are obtained, we may or may not be able to rule 
out certain values of /. 

One point which is worth emphasizing is the correlation of 7 with /. This study clearly shows 
that large values of / require smaller values of 7. The reason that this is important is as follows. 
The allowed range of 7 for a particular value of / is obtained from a fit to all CKM data, even 
those measurements which are unaffected by the presence of super symmetry. Now, the size of 7 
indirectly affects the branching ratio for B —>■ A s 7: a larger value of 7 corresponds to a smaller 
value of I Vt s | through CKM unitarity. And this branching ratio is among the experimental data 
used to bound SUSY parameters and calculate the allowed range of /. Therefore, the above 
7-/ correlation indirectly affects the allowed values of / in a particular SUSY model, and thus 
must be taken into account in studies which examine the range of /. For example, it is often 
the case that larger values of / are allowed for large values of 7. However, as we have seen 
above, the CKM fits disfavour such values of 7. 
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f 




h> 


/ 


(a B sy) 

yiX,/J, JJcent 


f = (SM) 
/ = 0.2 
/ = 0.4 
/ = 0.75 


65° - 123° 
70° - 129° 
75° - 134° 
86° - 141° 


16° - 35° 
16° - 34° 
15° - 34° 
13° - 33° 


36° - 97° 
32° - 90° 
28° - 85° 
23° - 73° 


(93°, 25°, 62°) 
(102°, 24°, 54°) 
(110°, 23°, 47°) 
(119°, 22°, 39°) 



Table 2: Allowed 95% C.L. ranges for the CP phases a, (3 and 7, as well as their central 
values, from the CKM fits in the SM (/ = 0) and supersymmetric theories, characterized by 
the parameter / defined in the text. 



/ 


sin 2a 


sin 2(3 


sin 2 7 


/ = (SM) 
/ = 0.2 
/ = 0.4 
/ = 0.75 


-0.91 - 0.77 
-0.98 - 0.65 
-1.00 - 0.50 
-1.00 - 0.14 


0.53 - 0.94 
0.52 - 0.93 
0.49 - 0.93 
0.45 - 0.91 


0.35 - 1.00 
0.28 - 1.00 
0.22 - 0.99 
0.16-0.91 



Table 3: Allowed 95% C.L. ranges for the CP asymmetries sin 2a, sin 2(3 and sin 2 7, from the 
CKM fits in the SM (/ = 0) and supersymmetric theories, characterized by the parameter / 
defined in the text. 



For completeness, in Table [3] we present the corresponding allowed ranges for the CP asym- 
metries sin 2a, sin 2(3 and sin 2 7. Again, we see that the allowed range of sin 2(3 is largely 
independent of the value of /. On the other hand, as / increases, the allowed values of sin 2a 
become increasingly negative, while those of sin 2 7 become smaller. 

The allowed (correlated) values of the CP angles for various values of / can be clearly seen 
in Figs. |4| and [|. As / increases from (SM) to 0.75, the change in the allowed sin2a-sin2/3 
(Fig. U) and a^j (Fig. [5|) regions is quite significant. 



In Sec. 2.1, we noted that (V^b/V^l, Bk and fB d y B Bd are very important in defining the 
allowed region in the p-r\ plane. At present, these three quantities have large errors, which 
are mostly theoretical in nature. Let us suppose that our theoretical understanding of these 
quantities improves, so that the errors are reduced by a factor of two, i.e. 









v cb 




B K 


fB d V 


'B Bri 



0.093 ±0.007 , 
0.94 ±0.07 , 

215 ± 20 MeV . (55) 

How would such an improvement affect the SUSY fits? 

We present the allowed 95% C.L. regions (/ = 0, 0.2, 0.4, 0.75) for this hypothetical situation 
in Fig. [?]. Not surprisingly, the regions are quite a bit smaller than in Fig. || More importantly 
for our purposes, the regions for the different values of / have become more separated from one 
another. That is, although there is still a region where all four / values are allowed, precise 
measurements of the CP angles have a better chance of ruling out certain values of /. 
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= 215 + 20MeV, B =0.94+0.07 
K 




Figure 7: Allowed 95% CL. region in p-rj space in the SM and in SUSY models, from a fit to 
the data given in Table [1], with the (hypothetical) modifications given in Eq. (|55|). From left 
to right, the allowed regions correspond to / = (SM, solid line), / = 0.2 (long dashed line), 
/ = 0.4 (short dashed line), / = 0.75 (dotted line). 

In Table |] we present the allowed ranges of a, (3 and 7, as well as their central values, for 
this scenario. Table |5| contains the corresponding allowed ranges for the CP asymmetries sin 2a, 
sin 2/? and sin 2 7. The allowed sin 2a-sin 2/5 and a-7 correlations are shown in Figs. [5| and |^, 
respectively. As is consistent with the smaller regions of Fig. ||, the allowed (correlated) regions 
are considerably reduced compared to Figs. [| and ||. As before, although the measurement of 
/3 will not distinguish among the various values of /, the measurement of a or 7 may. 

Indeed, the assumed reduction of errors in Eq. (|55"D increases the likelihood of this happening. 
For example, consider again Table 0, which uses the original data set of Table [TJ. Here we see 
that 65° < a < 123° for / = and 86° < a < 141° for / = 0.75. Thus, if experiment finds 
a in the range 86°-123°, one cannot distinguish the SM (/ = 0) from the SUSY model with 
/ = 0.75. However, consider now Table [|, obtained using data with reduced errors. Here, 
67° < a < 116° for / = and 97° < a < 137° for / = 0.75. Now, it is only if experiment finds 
a in the range 97°-116° that one cannot distinguish / = from / = 0.75. But this range is 
quite a bit smaller than that obtained using the original data. This shows how an improvement 
in the precision of the data can help not only in establishing the presence of new physics, but 
also in distinguishing among various models of new physics. 



4 Conclusions 

In the very near future, CP-violating asymmetries in B decays will be measured at 5-factories, 
HERA-B and hadron colliders. Such measurements will give us crucial information about the 
interior angles a, (3 and 7 of the unitarity triangle. If we are lucky, there will be an inconsistency 
in the independent measurements of the sides and angles of this triangle, thereby revealing the 
presence of new physics. 
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Figure 8: Allowed 95% C.L. region of the CP- violating quantities sin 2a and sin 2/3, from a fit 
to the data given in Table [1], with the (hypothetical) modifications given in Eq. (|55|). The upper 
left plot (/ = 0) corresponds to the SM, while the other plots (/ = 0.2, 0.4, 0.75) correspond 
to various SUSY models. 



/ 


a 




7 


(a, f3, 7) cent 


/ = (SM) 
/ = 0.2 
/ = 0.4 
/ = 0.75 


67° - 116° 
74° - 124° 
83° - 130° 
97° - 137° 


20° - 30° 
19° - 29° 
18° - 29° 
16° - 28° 


42° - 90° 
36° - 82° 
31° - 73° 
26° - 59° 


(93°, 24°, 63°) 
(102°, 24°, 54°) 
(110°, 23°, 47°) 
(119°, 22°, 39°) 



Table 4: Allowed 95% C.L. ranges for the CP phases a, (3 and 7, as well as their central 
values, from the CKM fits in the SM (/ = 0) and supersymmetric theories, characterized by 
the parameter / defined in the text. We use the data given in Table [1], with the (hypothetical) 
modifications given in Eq. (|55|). 
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Figure 9: Allowed 95% C.L. region of the CP-violating quantities a and 7, from a fit to the 
data given in Table [I], with the (hypothetical) modifications given in Eq. (|55|). The upper left 
plot (/ = 0) corresponds to the SM, while the other plots (/ = 0.2, 0.4, 0.75) correspond to 
various SUSY models. 



/ 


sin 2a 


sin 2(3 


sin 2 7 


/ = (SM) 
/ = 0.2 
/ = 0.4 
/ = 0.75 


-0.80 - 0.71 
-0.93 - 0.53 
-0.99 - 0.23 
-1.00 - -0.23 


0.64 - 0.86 
0.61 - 0.85 
0.57-0.85 
0.52 - 0.83 


0.44 - 1.00 
0.34 - 0.98 
0.27-0.91 
0.19 - 0.73 



Table 5: Allowed 95% C.L. ranges for the CP asymmetries sin 2a, sin 2(3 and sin 2 7, from the 
CKM fits in the SM (/ = 0) and supersymmetric theories, characterized by the parameter / 
defined in the text. We use the data given in Table [1], with the (hypothetical) modifications 
given in Eq. (|55|). 
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If present, this new physics will affect B decays principally through new contributions to 
B°-B° mixing. If these contributions come with new phases (relative to the SM), then the CP 
asymmetries can be enormously shifted from their SM values. In this case there can be huge 
discrepancies between measurements of the angles and the sides, so that the new physics will 
be easy to find. 

In fact, there are several models which have new contributions, with new phases, to B°-B° 
mixing. However, these models do not predict what those phases are. That is, there are no 
predictions for the values of CP asymmetries in such models. All that can be said is that large 
effects are possible. 

A more interesting possibility, from the point of view of making predictions, are models 
which contribute to B°-B° mixings and |e|, but without new phases. One type of new physics 
which does just this is supersymmetry (SUSY). There are some SUSY models which do contain 
new phases, but they suffer from the problem described above: lack of predictivity. How- 
ever, there is also a large class of SUSY models with no new phases. In this paper we have 
concentrated on these models. 

There has been an enormous amount of study of SUSY models over the past two decades. 
Much of this work has concentrated on the minimal supersymmetric standard model (MSSM), 
in which the new phases are constrained by limits on the electric dipole moments of the neu- 
tron and electron to be essentially zero. Taking into account supersymmetry breaking, in 
supergravity (SUGRA) models, the SUSY-breaking parameters of the MSSM are assumed to 
take a simple form at the Planck scale. However, there are a variety of ways to do this, so 
that in fact there is a fairly large class of SUSY models in which there are new contributions 
to B°-B° mixing, but no new phases. 

In this paper we have examined the predictions of such models for the unitarity triangle and 
explored the extent to which this type of new physics can be discovered through measurements 
of the sides and angles of the unitarity triangle. 

In these models, there are new, supersymmetric contributions to K°-K°, B®-B® and B®- 
B® mixing. The key ingredient in our analysis is the fact that these contributions, which add 
constructively to the SM, depend on the SUSY parameters in essentially the same way. That 
is, so far as an analysis of the unitarity triangle is concerned, there is a single parameter, /, 
which characterizes the various SUSY models within this class of models (/ = corresponds to 
the SM). For example, the values / = 0.2, 0.4 and 0.75 are found in minimal SUGRA models, 
non-minimal SUGRA models, and non-SUGRA models with EDM constraints, respectively. 

We have therefore updated the profile of the unitarity triangle in both the SM and some 
variants of the MSSM. We have used the latest experimental data on \V c b\, \V u b/V c b\, AMa 



£ s = fs d y BbJ fB s \j Bb s - In addition to / = (SM), we considered the three SUSY values of 
/: 0.2, 0.4 and 0.75. 

We first considered the profile of the unitarity triangle in the SM, shown in Fig. |^. At 
present, the allowed ranges for the CP angles at 95% C.L. are 



and AM S , as well as the latest theoretical estimates 





65° < a < 123 



O 



16° < (3 < 35 



O 



36° < 7 < 97° , 



(56) 
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or equivalently, 

- 0.91 < sin 2a < 0.77 , 0.52 < sin 2(3 < 0.94 , 0.35 < sin 2 7 < 1.00 . (57) 



We then compared the SM with the different SUSY models. The result can be seen in 
Fig. [| As / increases, the allowed region moves slightly down and to the right in the p-77 
plane. The main conclusion from this analysis is that the measurement of the CP angle (3 
will not distinguish among the SM and the various SUSY models - the allowed region of (3 is 
virtually the same in all these models. On the other hand, the allowed ranges of a and 7 do 
depend on the choice of /. For example, larger values of / tend to favour smaller values of 7. 
Thus, with measurements of 7 or a, we may be able to rule out certain values of / (including 
the SM, / = 0). However, we also note that there is no guarantee of this happening - at present 
there is still a significant region of overlap among all four models. 

Finally, we also considered a hypothetical future data set in which the errors on |V™&/14&|, 
Bx and fg d yBs d , which are mainly theoretical, are reduced by a factor of two. For two of 

these quantities (IV^/V^I and fs d \J B Bd ), this has the effect of reducing the uncertainty on the 
sides of the unitarity triangle by the same factor. The comparison of the SM and SUSY models 
is shown in Fig. [7|. As expected, the allowed regions for all models are quite a bit smaller than 
before. Furthermore, the regions for different values of / have become more separated, so that 
precise measurements of the CP angles have a better chance of ruling out certain values of /. 



Acknowledgements: 



We would like to especially thank Seungwon Baek, John Ellis, Toru Goto and Pyungwon 
Ko for very helpful discussions, and Laksana Tri Handoko for numerical help. Thanks are 
also due to Hans-Gunther Moser, Olivier Schneider and Fabrizio Parodi for providing us the 
experimental information on the LEP measurements of the B-B mixings and instructions in 
implementing the resulting constraints. D.L. is grateful to G. Azuelos and I. Trigger for their 
invaluable help with MINUIT and PAW. This research was financially supported by NSERC 
of Canada and FCAR du Quebec. 

Note added: As we were submitting this paper for publication, we noted the preprint by F. 
Parodi, P. Roudeau and A. Stocchi (LAL 99-03, DELPHI 99-27 CONF 226, |hep-ex/9903063l) , 



in which the constraints on the parameters of the CKM matrix are presented in the SM. Though 
our input parameters differ, in particular the assumed errors, our results in the SM are similar 
to theirs but not identical. 



References 

[1] N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531; M. Kobayashi and K. Maskawa, Prog. Theor. 
Phys. 49 (1973) 652. 

[2] L. Wolfenstein, Phys. Rev. Lett. 51 (1983) 1945. 



27 



[3] For reviews, see, for example, Y. Nir and H. R. Quinn in B Decays, edited by S. Stone 
(World Scientific, Singapore, 1994), p. 362; I. Dunietz, ibid., p. 393; M. Gronau, Proceed- 
ings of Neutrino 94, XVI International Conference on Neutrino Physics and Astrophysics, 
Eilat, Israel, May 29 - June 3, 1994, eds. A. Dar, G. Eilam and M. Gronau, Nucl. Phys. 
(Proc. Suppl.) B38 (1995) 136. 



[4] 
[5] 

[6] 
[7] 
[8 
[9 

[io; 
[ii 



CO. Dib, D. London and Y. Nir, Int. J. Mod. Phys. A6 (1991) 1253. 

For a review of new-physics effects in CP asymmetries in the B system, see M. Gronau 
and D. London, Phys. Rev. D55 (1997) 2845, and references therein. 

M. Ciuchini, G. Degrassi, P. Gambino and G.F. Giudice, Nucl. Phys. B534 (1998) 3. 

F. Krauss and G. Soff, preprint |hep-ph/9807238| . 

C. Caso et al. (Particle Data Group), Eur. Phys. J. C3 (1998) 1. 

M. Neubert, in Heavy Flavors, Second Edition, edited by A.J. Buras and M. Lindner 
(World Scientific, Singapore, 1997) and |hep-ph/9702375| . 

A. Czarnecki and K. Melnikov, Phys. Rev. Lett. 78 (1997) 3630; Nucl. Phys. B505 (1997) 
65. 



M. Shifman, in Proc. PASCHOS 1995, Baltimore, Maryland, 22-25 March 1995 and |hep^ 
ph/9505289 



[12] F. Parodi, in Proc. of XXIXth. Int. Conf. on High Energy Physics, Vancouver, B.C., 1998. 
[13] S. Mele, preprint CERN-EP/98-133, |hep-ph/9810333| (1998). 

[14] A.J. Buras, W. Slominski, and H. Steger, Nucl. Phys. B238 (1984) 529; ibid. B245 (1984) 
369. 

[15] T. Inami and C.S. Lim, Progr. Theor. Phys. 65 (1981) 297. 

[16] S. Herrlich and U. Nierste, Nucl. Phys. B419 (1994) 292. 

[17] A.J. Buras, M. Jamin and PH. Weisz, Nucl. Phys. B347 (1990) 491. 

[18] S. Herrlich and U. Nierste, Phys. Rev. D52 (1995) 6505. 

[19] T. Draper, preprint |hep-lat/981006l (1998). 

[20] S. Sharpe, preprint |hep-lat/9811006| (1998). 

[21] J. Alexander, in Proc. of XXIXth. Int. Conf. on High Energy Physics, Vancouver, B.C., 
1998. 

[22] A. Ali and D. London, Z. Phys. C65 (1995) 431; Nucl. Phys. Proc. Suppl. 54A (1997) 
297. 

[23] C. Bernard et al, Phys. Rev. Lett. 81 (1998) 4812. 

28 



[24 
[25 

[26 
[27; 
[28 

[29 
[30 

[31 
[32 

[33 

[34 
[35 



[37; 
[38 



[40 
[41 
[42 



H.G. Moser and A. Roussarie, Nucl. Instr. Meth. A384 (1997) 491. 

Seminar presented by P. Shawhan (KTEV Collaboration), Fermilab, IL, Feb. 24, 1999; see 
also |http:/ /fnphyx-www.fnal.gov/experiments/ktev/epsprime/| . 

G.D. Barr et al. (NA31 Collaboration), Phys. Lett. B317 (1993) 233. 

L. Wolfenstein, Phys. Rev. Lett. 13 (1964) 380. 

A.J. Buras, M. Jamin and M.E. Lautenbacher, Nucl. Phys. B370 (1992) 69; ibid. B400 
(1993) 37; ibid. B400 (1993) 75; ibid. B408 (1993) 209. 

M. Ciuchini et al, Phys. Lett. B301 (1993) 263; Nucl. Phys. B415 (1994) 403; 

S. Bertolini, M. Fabbrichesi and J.O. Egg, Nucl. Phys. B499 (1995) 197; ibid. B476 (1996) 
225. 



A.J. Buras and L. Silvestrini, preprint TUM-HEP-334/98, [hep-ph/ 98 11471 



For recent reviews, see R. Gupta, preprint |hep-ph/9801412| and L. Conti et al., Phys. Lett. 
B421 (1998) 273. 

W.A. Bardeen, A.J. Buras and J.M. Gerard, Phys. Lett. B180 (1986) 133; Nucl. Phys. 
B293 (1987) 787; Phys. Lett. B192 (1987) 138; CO. Kdhler et al, Phys. Rev. D58 (1998) 
014017. 



A.J. Buras et al, preprint TUM-HEP-316-98, |hep-ph/980647T 



Y.-Y. Keum, U. Nierste and A.I. Sanda, preprint KEK-TH-612, FERMILAB-Pub-99/035- 
T, |hep-ph/990323q . 

A. Masiero and H. Murayama, preprint |hep-ph / 9903363 . 
CDF Collaboration, CDF/PUB/BOTTOM/CDF/4855. 
F. Parodi, P. Roudeau and A. Stocchi, preprint LAL 98-49, |hep-ph/9802~289 



M. Gronau and D. Wyler, Phys. Lett. B265 (1991) 172. See also M. Gronau and D. 
London, Phys. Lett. B253 (1991) 483; I. Dunietz, Phys. Lett. B270 (1991) 75. Improve- 
ments to this method have recently been discussed by D. Atwood, I. Dunietz and A. Soni, 
Phys. Rev. Lett. 78 (1997) 3257. 

R. Aleksan, I. Dunietz and B. Kayser, Z. Phys. C54 (1992) 653. 
See, for example, H.P Nilles, Phys. Rep. 110 (1984) 1. 

M. Dugan, B. Grinstein and L.J. Hall, Nucl. Phys. B255 (1985) 413; S. Dimopoulos and 
S. Thomas, Nucl. Phys. B465 (1996) 23. 



[43] A.G. Cohen, D.B. Kaplan and A.E. Nelson, Phys. Lett. B388 (1996) 588; A.G. Cohen, 
D.B. Kaplan, F. Lepeintre and A.E. Nelson, Phys. Rev. Lett. 78 (1997) 2300. 

29 



[44] T. Falk, K.A. Olive and M. Srednicki, Phys. Lett. B354 (1995) 99; T. Falk and K.A. Olive, 
ibid. B375 (1996) 196. 



[45 
[46 

[47 
[48 



T. Nihei, Prog. Theor. Phys. 98 (1997) 1157. 

G.C. Branco, G.C. Clio, Y. Kizukuri and N. Oshimo, Phys. Lett. B337 (1994) 316; 
Nucl. Phys. B449 (1995) 483. 

R. Contino and I. Scimemi, preprint ROME1-1216-98, frep-ph/ 98091571 



D.A. Demir, A. Masiero and O. Vives, preprint SISSA/EP/140/98, IC/98/228, |hep- 



ph/9812337 . 



[49] S. Baek and P. Ko, preprint KAIST-20/98, SNUTP 98-139, [hep-ph/98 12229 . 

[50] D. Chang, W.-Y. Keung and A. Pilaftsis, Phys. Rev. Lett. 82 (1999) 900. 

[51] J. P. Silva and L. Wolfenstein, Phys. Rev. D55 (1997) 5331. 

[52] L. Randall and S. Su, Nucl. Phys. B540 (1999) 37. 

[53] M. Ciuchini et al, JHEP 9810: 008, 1998. 

[54] T. Goto, T. Nihei and Y. Okada, Phys. Rev. D53 (1996) 5233. 

[55] T. Goto, Y. Okada, Y. Shimizu and M. Tanaka, Phys. Rev. D55 (1997) 4273. 

[56] T. Goto, Y. Okada and Y. Shimizu, Phys. Rev. D58: 094006 (1998); KEK-TH-611 (1999) 
(in preparation). 

[57] T. Goto et al, preprint KEK-TH-608, KEK Preprint 98-206, |hep-ph/9812369| . 

[58] S. Bertolini, F. Borzumati, A. Masiero and G. Ridolphi, Nucl. Phys. B353 (1991) 591. 

[59] M. Dine and A.E. Nelson, Phys. Rev. D48 (1993) 1277; M. Dine, A.E. Nelson and 
Y. Shirman, Phys. Rev. D51 (1995) 1362; M. Dine, A.E. Nelson, Y. Nir and Y. Shir- 
man, Phys. Rev. D53 (1996) 2658. 

[60] M. Ciuchini, G. Degrassi, P. Gambino and G.F. Giudice, Nucl. Phys. B527 (1998) 21; 
P. Ciafaloni, A. Romanino and A. Strumia, Nucl. Phys. B524 (1998) 361; F. Borzumati 
and C. Greub, Phys. Rev. D58 (1998) 74004, |hep-ph/9809438| ; T.M. Aliev and E.O. Iltan, 
Phys. Rev. D58 (1998) 95014. 



30 



